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Abstract 



We continue our analysis of the general A^-Higgs-doublet model and focus of the Higgs 
potential description in the space of gauge orbits. We develop a geometric technique 
that allows us to study the global minimum of the potential without explicitly finding 
its position. We discuss symmetry patterns of the NHDM potential, and illustrate the 
general discussion with various specific variants of the three-Higgs-doublet model. 



> 

1 Introduction 



! The electroweak theory relies on the Higgs mechanism of the electroweak symmetry breaking 
(EWSB). Many variants of the Higgs mechanism have already been suggested, see e.g. pp. 
They range from conceptually simple ones where one just adds extra particles to the scalar 
sector of the Standard model to the models involving new gauge interactions, extra dimensions 
at the TeV scale, supersymmetry, etc. It is not presently known what particular variant is 
realized in Nature, and it is the key goal of the LHC and of future linear colliders to answer 
that question. However, in order to safely interpret possible manifestations of the physics 
beyond the Standard model (BSM) in the future experimental data, it is necessary to clearly 
understand in advance what each model is capable of. 

There is a broad class of BSM models known as A^-Higgs-doublet models (NHDM), in 
which one gets New Physics just by introducing several copies of a Higgs doublet with the 
usual electroweak quantum numbers. They are among the most conservative extensions of the 
Standard model and are driven by the idea to take the notion of "generations" to the scalar 
sector. This idea is quite natural because, first, there is no strong argument that would favor 
a single Higgs doublet, second, we do have generations of fermions without understanding 
their origin, and third, since fermion masses and mixing patterns result from their Yukawa 
interactions with scalars, there might exist a link between generations of scalars and fermions. 
Additional support is lent by supersymmetric models and certain low-energy realizations of 
superstring/brane models, which require presence of several Higgs generations, see [2] and 
references therein. 
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Multi-Higgs-doublet models were studied a lot in the past decades (see some references 
below). Virtually all the analyses were restricted to simple variants of the Higgs self- interaction 
potential, with interesting phenomenological consequences observed. However one cannot 
know if these particular cases exhaust all the New Physics offered by several doublets or if 
there are some "hidden treasures" lurking in the multi-dimensional parameter space of the 
model. To find them, one certainly needs to study the iV-doublet models with most general 
scalar potential and Yukawa sector. 

Algebraically, working with the most general iV-Higgs-doublet potential is a hard task. 
Difficulties arise even at the tree level within the scalar sector of the theory. A very repre- 
sentative case is given by the two-Higgs-doublet model (2HDM), [U EJ H]. Here, the Higgs 
potential cannot be minimized explicitly in the general case. As a result, for a long time only 
relatively simple variants of the 2HDM remained analyzed, while the most general 2HDM 
stood barely touched. However, in the last several years new tools have been developed which 
led to many insights into the properties of the general 2HDM. These methods are based on 
the idea of the reparametrization invariance of the physical observables in the scalar sector: a 
unitary transformation between the Higgs doublets changes the parameters of the lagrangian, 
but the physical observables such as the Higgs boson masses and the sum of squares of the 
vacuum expectation values remain the same. This idea was implemented via the tensorial for- 
malism at the level of Higgs fields [5j EJ [7J EJ [9] or via geometric constructions in the space of 
gauge- invariant bilinears [TU 1 fill fT2l [T3l 11%] . In the latter case the formalism was extended to 
include non- unitary reparametrization transformations [151 EE] EE] , which revealed interesting 
geometric properties of the 2HDM in the orbit space equipped with the Minkowskian metric. 

This formalism led to discovery of new features of the scalar sector of 2HDM that could 
not be found or were just too obscure to be observed with more tranditional techniques. They 
include, among others, the following results: 

• There can exist at most two distinct minima of the Higgs potential, [T6] : 




• charge-breaking and neutral minima never coexist (although this fact was known since 
[18] . its proof becomes straightforward in the Minkowskian formalism, [15]); 

• CP- conserving and CP- violating minima cannot coexist, [15J; 

• there are exactly six conjugacy classes of symmetries that can be imposed on the scalar 
lagrangian [T21 [T3J [Hj, all specific realizations of each symmetry can be mapped onto 
each other by a reparametrization transformation. Remarkably, this conclusion holds 
even for the scalar potential alone, despite the fact that it can be more symmetric than 
the full scalar lagrangian, [16J; 

• the maximal breaking of a discrete symmetry of the potential consists in removing just 
a single Z 2 factor, [16J; 

• the full tree-level phase diagram of the scalar sector was described and possible phase 
transitions were identified [16]; 

• the thermal evolution of this phase diagram was calculated in the first non-trivial ap- 
proximation, with interesting results for possible multiple phase transition paths of the 
early Universe, [19]. 
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Apart from that, many algebraic conditions were derived in a basis-invariant form (positivity 
constraints, conditions for existence of any specific symmetry in a given lagrangian, conditions 
when this symmetry is broken). 

It is a natural idea to extend these successful techniques to A" doublets. The general 
NHDM is obviously more involved than 2HDM, both at the level of scalar sector and Yukawa 
interactions, see examples in (20j [2TJ [22] . Some properties of the general NHDM were analyzed 
in pH 123 HU 125J |2Sl [21], with a special emphasis on CP- violation, [Ml [28] . However, a 
method to systematically explore all the possibilities offered with N doublets was still missing. 

It is the aim of this paper to fill this gap, at least partially. Making use of the results of the 
companion paper [29] , we study here the general NHDM Higgs potential in the space of bilin- 
ears. We develop a geometric technique that allows us to study the properties of the minima 
of the Higgs potential without explicitly finding its position. The general constructions are 
illustrated with several variants of the three-Higgs-doublet model (3HDM) displaying various 
patterns of spontaneous symmetry breaking. Note that in this paper we focus only on the 
Higgs potential; the interaction between the Higgs sector and the fermions is an intriguing 
topic, which we do not touch here. 

The paper is organized as follows. In Section 2 we briefly overview some of the results 
obtained in [29J concerning the orbit space of the NHDM. Then, in Section 3 we describe 
the equipotential surface technique. Section 4 is devoted to symmetries of NHDM. These 
general discussions are used in Section 5 to construct several variants of 3HDM with various 
symmetries and symmetry breaking patterns. We conclude the paper with a discussion of the 
results and an outline for future work. 



A detailed description of the orbit space of NHDM was given in [29J. Here we just introduce 

the notation and briefly summarize some of its results. 

In the Af-Higgs-doublet model we introduce A" Higgs doublets a , a = 1, . . . , N . The gen- 
eral renormalizable Higgs potential of NHDM is constructed from the gauge-invariant combina- 
tions which describe the gauge orbits in the Higgs space. Following [TTJ [121, H3l I2S] , we 
organize them in the i^-matrix: K a ^ = (0^0 a ), which is a hermitian and positive-semidefinite 
N x N matrix with rank < 2. A rank-2 i^-matrix corresponds to a charge-breaking minimum, 
while a rank-1 matrix corresponds to a neutral minimum. Any .ff -matrix can be brought to a 
diagonal form by an appropriate reparametrization transformation (a unitary or antiunitary 
rotation among the doublets). 

The .fT-matrix can be decomposed as 



Here ljy is the N x N unit matrix and Aj are generators of the algebra su(N). The coefficient 
in front of the unit matrix in (Cp) is chosen for future convenience. Eq. (CQ) defines a linear and 



2 The orbit space of NHDM 



2.1 General N 
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invertible map from the space of ^-matrices to the space of real vectors^ = (r ,Tj) G R^ 2 , 
which we call the adjoint space. As it was noted in [121 122], this map is not surjective. 
Therefore, the orbit space does not cover the entire space of vectors r M , but is represented by 
a certain manifold embedded in it, which we denote by V$. This manifold can be described 
algebraically with a list of polynomial (in) equalities in components of r M . This list starts with 

N — 2 

r >0, r 2 -r*>0, d ijk r irj r k + —===== (r 2 - 3#*) = , (2) 

where dijk are the fully symmetric SU(N) tensors, see details in [29]. The manifold of the 
neutral vacua can be obtained from (j2J) by turning the second inequality into equality. In 
2HDM we have just the first two inequalities, which define the surface and interior of the 
forward lightcone in R 4 . In 3HDM the third equality is added, which carves out a specific 
8-dimensional manifold in the forward lightcone in R 9 . 

Some geometric properties of the orbit space manifold V$ were established in |29j. The 
dimensionalities of the charge-breaking and neutral parts of the orbit space manifold, Vc and 
Vn, are 4N — 4 and 2N — 1, respectively. Since all the equations in (J2J) are homogeneous, the 
orbit space has a conical shape in the adjoint space. Its base (i.e. a r$ = const section) can 
be described in the (N 2 — l)-dimensional vector space of n = f/ro, where Vn becomes the 
complex projective space CP^ -1 embedded in R^ _1 , while the entire orbit space V$ can be 
reconstructed as the symmetric join of Vat (the union of all points that lie on line segments 
drawn between all possible pairs of Vat points). 

It is also very useful to consider the root space of NHDM. It is the (N — l)-dimensional 
subspace spanned by the Cartan subalgebra of su(N), i.e. the space of coordinates rij in 
front of the diagonal matrices Aj in (pQ), with all the other coordinates set to zero. Any point 
of the orbit space can be brought to the root space by an appropriate reparametrization 
transformation. In the root space, the orbit space is represented by an (N — l)-simplex: the 
neutral orbit space corresponds to the N vertices of an N- dimensional regular tetrahedron, 
while the charge-breaking orbit space corresponds to the iV(iV — l)/2 edges joining all pairs 
of vertices. Note that the interior of the tetrahedron does not belong to the orbit space, so 
that a curious feature of the NHDM is the existence of a "hole" inside the orbit space (or an 
empty "inner cone" in the r M space). 



2.2 The case of 3HDM 

Let us consider more closely the orbit space of the three-Higgs-doublet model (3HDM). For 
N = 3, the matrices Aj, i = 1, . . . , 8, are the standard Gell-Mann matrices, and the Cartan 
subalgebra is formed by the diagonal matrices 

A 3 = diag(l, — 1, 0) and A 8 = — ^=diag(l, 1, — 2). (3) 

v3 



1 The notation r p alludes to the Minkowski-space formalism of [T5J [TBJ [T31 [TT] . Since we limit ourselves only 
to the (anti)unitary reparametrization transformations, we will not use this formalism in this paper; is just 
a short notation for (ro,Ti). 
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The explicit expressions for the field bilinears are: 
1 



r n = —=TrK 



V3 



(4) 



r 3 



402) 



r 8 



+ 



-2( 



2^ 



ri = Re((^0 2 ) , r 2 = Im(</>j</> 2 ) , r 4 = Re(</40 3 ) , 
r 5 = lm(4M , r 6 = Re(0 2 3 ) , r 7 = lm(0 2 3 ) . 

It is also useful to group the last six real coordinates into three "complex coordinates" : 

ru = n + ir 2 , r 45 = r 4 - ir b , r 67 = r 6 + ir 7 , (5) 

with the same labels for the normalized vectors 7ij = Tj/r . 

The root space of 3HDM is given by the (n 3 , n§) plane, shown in Fig. [TJ On the root plane, 
the orbit space is represented by three points P, P', P" (neutral manifold) and the three line 
segments joining them (the charge-breaking manifold). 




Figure 1: The root plane for 3HDM. Shown are the unit circle (section of the lightcone), the 
three points P, P', P" representing the neutral manifold, the three line segments representing 
the charge-breaking manifold, as well as the three sets of coordinates: (n 3 ,ng) shown in black 
solid lines, (n' 3 ,n' 8 ) shown in dashed blue lines, and (n 3 ,n 8 ) shown in dash-dotted red lines. 

The full orbit space V$ is symmetric under the group S3 of permutations of the three 
doublets, but the description in terms of rii breaks this symmetry. In order to restore it, we 
introduce two extra sets of the coordinates, n! i and n", which are obtained from rti by cyclic 
permutations of the three doublets 4>i — > 02 — > 03 — > 0i an d its inverse: 



Tin 



1 Vs 



Tin = — 



- 2 n, - —n 8 



ru 



-n 3 



12 



n 67 , n 



45 



"12 , n', 



67 



V3 1 „ „ 

—n 3 - -n 8 , n 12 = n 45 , n 45 



"67, n 



67 



n 45 



n 12 . (6) 



5 



Note that n 3 + n' 3 + n 3 = rig + n' s + n'l = 0. In terms of these coordinates, the three neutral 
points are 



P : n 3 = 0, n 8 = -1 , P' : n' 3 = 0, n' 8 = -1 
Finally, we also introduce the baricentric coordinates 



P" : n'i = 0, n' 8 ' 



-1. 



1 - 2n« 



P 



, 1 " 2n' 8 



1 - 2n" 



with p + p' + p" = 1 , 



(7) 



(8) 



which are proportional to the distances from a given point to the three edges of the triangle. 
Using all these coordinates, one can write the i^-matrix of 3HDM in a very symmetric way: 



K a b = r 



^12 ^19 ^ 



n 12 V% 



" n'* 



\ 



n 



12 



n 



(9) 
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Let us now turn to the third algebraic conditions in ()2]), which we will refer to as "the 
rf-condition" : 

r— 3n 2 - 1 

V3d ijk ninjn k = — - — . (10) 

With the three sets of coordinates, this condition can be rewritten symmetrically as 



p\ni 2 \ +p\n 12 \ +p \n 12 \ -3pp'p" -=Re(n 12 n 12 n" 2 ) = . 

v3 

For the neutral manifold, the <i-condition can be simplified further. In this case we have 



n I II I 1 2 o I I 

6pp = \ni2\ , opp = \n 



11 1 2 

12 I 7 



Sp p = \n 



I |2 
12 1 



(12) 



Denoting the sum of the phases of n 12 , n' 12 , and n" 2 by 7, one can cast the ci-condition for the 
neutral orbit space into 

pp'p"{l -cos 7) = 0. (13) 



3 Minimization of the general NHDM Higgs potential 

The key attractive feature of the adjoint space formalism is the simplicity with which one can 
describe the Higgs potential and discuss its symmetries. The most general gauge-invariant 
and renormalizable Higgs potential with iV Higgs doublets can be written in a very compact 
way 

V = -M^ + X -k^r v . (14) 

As we mentioned above, we do not use the Minkowski-space formalism of [151 [TB| IT3"t [T7] in 
this paper, so r M just stands for (r , r^), % — 1, . . . , iV 2 — 1, and all index contractions are meant 
to be Euclidean: 

1 1 

V = -M r - Min + -A 00 rl + A^ron + ^r.r^ . (15) 
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Note that when constructing a potential, we have the full freedom to choose iV 2 components 
of M M and N 2 (N 2 + l)/2 components of A^. 

The first thing one needs to do with the potential is to minimize it. Unfortunately, this 
task cannot be done with straightforward algebra for a general potential. Alternatively, one 
can start from a minimum and construct a potential, which by construction has a minimum 
at a chosen point. In this case, however, one must be sure that the potential does not possess 
another, deeper minimum, which again cannot be verified with straightforward algebra. 

In [16] a geometric technique was suggested for 2HDM that allowed one to study the 
properties of the global minima of the Higgs potential without knowing its exact position. 
With this technique, one constructs equipotential surfaces and finds their intersection with 
the orbit space. Although this approach did not help to explicitly minimize a given potential, 
it offered a good understanding of the entire spectrum of possibilities (symmetries and their 
breaking, coexistence of minima, phase transitions), which can in principle arise at the tree- 
level in 2HDM. 

In this Section we generalize this approach to A" doublets. 

3.1 Equipotential surfaces 

The Higgs potential (JT3J) is a quadratic form in the components of r M . In order to be physically 
realizable in terms of doublets, r M must lie in the orbit space, whose complicated shape we 
briefly described in Section Let us, however, put aside this requirement for a moment and 
consider the quadratic form f JT5|) in the entire adjoint space M. n2 . 

We choose a real number C and construct an equipotential surface Aic as the set of all 
points satisfying V(r^) = C . Taking another value of C will result in another equipotential 
surface constructed from the same potential, and so on. Several observations can be inferred 
from this definition. 

• Each Aic is an (A^ 2 — l)-quadric, i.e. a second-order algebraic manifold embedded 
in R N : a hyperboloid, an ellipsoid, etc. If the potential has flat directions, then the 
equipotential surfaces are cylindrical along these directions. 

• Two equipotential surfaces of the same potential but with different Cj do not intersect. 

• The equipotential surfaces are linearly ordered: if C\ < C 2 < C 3 , then Aic 2 hes between 
Aid an d M-c 3 - One can also say that Aic 2 ^ es higher than Aid, but lower than Aic 3 - 

Thus, a given Higgs potential can be associated with a nested family of equipotential surfaces, 
which covers the entire adjoint space. There are two sorts of families of equipotential surfaces 
to be considered: 

• .R-family: solutions to equation V(r) = C exist for any real C. 

• i? + -family: solutions to equation V(r) = C exist only for all real C larger or equal to 
some value. 

Which case we have depends solely on the signs of the Euclidean eigenvalues of A MI/ . The 
quadratic form 

A^rV = A 00 r 2 + 2A 0i r r i + A, (16) 
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is real and symmetric, therefore, it can be diagonalized in the (r , r) space by some rotation of 
the coordinates: (r , r) — > {r q }, q = 1, . . . , iV 2 . After this transformation, the quadratic form 
f|T6|) becomes A q r 2 , with at least one positive eigenvalue among A q (otherwise, the potential 
would not be bounded from below in the orbit space). The i? + -family corresponds to the case 
when all A q > 0, that is when A M „ is positive semidefmite in the entire adjoint space M N . The 
i?-family corresponds to situations when at least one of A q is negative. 

Note that we do not claim that a rotation diagonalizing A^ u is physical and can be realized 
by a certain transformation of doublets. It well may be not because there is no transforma- 
tion of the doublets that mixes r and keeping r 2 . + rf constant. We use this unphysical 
transformation just to characterize the spectrum of the matrix A^ u and define the shape of 
the equipotential surfaces. 

3.2 Positivity conditions 

The physical Higgs potentials we consider must be stable, i.e. bounded from below. Potentials 
with a quadratic V2 and a quartic V4 terms can be stable in two senses. Using terminology 
of [12], we say that a potential is stable in a strong sense if its V4 increases along any ray in 
the space of Higgs fields. A potential stable in a weak sense is the one whose V4 can have flat 
directions, but Vj> grows along them. In this paper we focus on potentials stable in a strong 
sense. Potentials stable in a weak sense can be analyzed in a similar manner. 

An i? + -family of equipotential surfaces automatically guarantees the stability of the po- 
tential. For an i?-family, we require that the quartic part of the potential be positive definite 
everywhere in the orbit space (except the origin): 

A^r'V^ > for all non-zero in the orbit space . (17) 

Unfortunately, we could not find a compact algebraic reformulation of this criterion in terms 
of components of A^ only (an attempts to find such conditions using hyperspherical coordi- 
nates can be found in [27]). It is possible that such a criterion will involve its Minkowskian 
eigenvalues, similarly to the result [T5| [T7] for 2HDM, although a sufficient condition is easy 
to formulate for NHDM as well, see [25] . 

On the other hand, the positivity constraint on the Higgs potential can be easily described 
geometrically, which helps visualize the freedom we have when constructing NHDM potentials. 
Consider again the quadratic form A^ lu r fJ- r l/ in the entire M. N . Clearly, it stays sign-definite 
and just increases in absolute value along any ray drawn from the origin. Consider now the 
part of R N2 filled by rays where A^r^r" < 0. The potential is stable if and only if this 
"negative region" does not intersect the orbit space. 

This description is illustrated schematically in Fig. [2] on a generic (ro, rj)-plane. Here, the 
neutral orbit space is indicated by the thick rays, the charge-breaking orbit space is shown by 
the gray region, while the "negative region" is the dashed area labeled as A < 0. An interesting 
possibility absent in 2HDM is that the "negative region" can be placed not only outside, but 
also inside the lightcone, see Fig. |2j right. In this case the negative region is contained in the 
empty inner cone. An example of 3HDM potential which realizes such an exotic situation is 

V = -M r - M iri + ^A 00 r 2 + -A^r^ , (18) 
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orbit space 




Figure 2: Location of the negative region (shown as dashed area and labeled as A < 0) with 
respect to the orbit space on a generic (ro,rj) plane. The neutral orbit space is shown with 
thick lines, the charge-breaking orbit space is shown as a gray area. For 3HDM, the negative 
region can lie either outside the lightcone (left) or inside the empty inner cone (right). 



with a positive-definite and a negative A 00 satisfying 

n . K,,min 
U > A 00 > , 

where Aj jmin is the smallest eigenvalue of Ay. Note also that it is a specific property of 3HDM 
that the negative region can be placed either inside or outside the orbit space, because the 
orbit space of 3HDM isolates the inner cone from the region outside the lightcone. In NHDM 
with N > 3 the orbit space is much more "perforated", so that the inner cone is in fact 
connected with the outer region. Thus, intermediate situations between the two extremes 
shown in Fig. [2] are also possible. 



3.3 Searching for the global minimum 

With the help of equipotential surfaces, the search for the global minimum of the potential 
can be reformulated in a very intuitive geometric way. 

Suppose that we have made sure that the positivity conditions are satisfied. Then, an 
equipotential surface M. c with a sufficiently negative C does not intersect the orbit space. 
Let us now increase the value of C. The lowest value of C = Co, at which we have the first 
contact with the orbit space, defines the critical equipotential surface Aic - The contact point 
or points give the position(s) of the global minimum of the potential. Thus, in order to identify 
the global minimum of a given potential, we need to find the unique equipotential surface that 
only touches but does not intersect the orbit space. This fact, as well as the uniqueness of the 
critical equipotential surface, follow from the fact that the orbit space is a connected manifold 
and from the nesting property of the equipotential surfaces. 

This construction is illustrated in Fig. Son a generic plane (r ,rj). The outer boundary 
of the orbit space is shown schematically as a pair of rays. Four equipotenial surfaces M.\ 
through M.^ of a Higgs potential are plotted as pairs of hyperboles. The equipotential surface 
M.\ corresponds to such a low value of C that it never intersects the orbit space, so it is not 
realized in terms of doublets. M. 2 (shown in bold) is the lowest equipotential surface that 
touches the orbit space, its contact point A giving the position of the global minimum. Any 
equipotential surface lying higher, such as and .M 4 , necessarily intersect the orbit space 
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Figure 3: Search for the global minimum via equipotential surfaces, see text for explanation. 

along some region, not just a single point. However sometimes it can also lead to additional 
local minima, such as point B on Ai^. 

Although this geometric criterion for identification of the global minimum does not help 
to minimize a given potential, it is of much use when we construct potentials with a specific 
explicit symmetry and a specific pattern of spontaneous symmetry breaking. It will guarantees 
us that we are indeed working with the global minima of the potential. 

4 Symmetries of NHDM 

One of the main motivations to study non-minimal Higgs sectors is a possibility to introduce 
new symmetries via the scalar sector. In particular, three questions are of much importance: 

1. Which symmetries can the lagrangian of a given non-minimal Higgs sector, e.g. NHDM, 
have? 

2. What are the possible ways to spontaneously break these symmetries? Which phenom- 
ena accompany symmetry breaking? 

3. How can these symmetries be extended into the fermion sector? 

A thoroughly studied example is again given by the 2HDM. Various symmetries have been 
discussed starting from the original paper on the subject [3], for a recent review see [I]. Since 
various symmetries with the same group could be mapped to each other by an appropriate 
reparametrization transformation, they represent the same symmetry of the Higgs lagrangian 
and might differ only in the fermionic sector. Therefore, if one focuses on the scalar lagrangian 
only, one should look at conjugacy classes of symmetries, not at their particular implementa- 
tions, see [26] for a discussion of this point. 

In [16] the complete classifications of symmetry classes in the most general 2HDM was 
given. Only six conjugacy classes of symmetry can exist in 2HDM, with groups Z2, {Z2) 2 , 
(Z2) 3 , Z2 x 0(2), and 0(3). Reparametrization invariant criteria of the presence of a (hidden) 
symmetry were given in [HI H2J [131 El EEl EE E], with a particular emphasis on the CP- 
violation. Which of these symmetries can be extended into the fermionic sector has also 
received some attention recently, see [301 EB [32] . 
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In NHDM, one can expect an even richer spectrum of possible symmetries and symmetry 
breaking; some of them were studied inj2Qll2Ill22ll2Hll23llHII31l2Sll2S]- Nevertheless, despite 
all these efforts, we are still far from being able to write down the list of all symmetries possible 
in NHDM. 

In the Section we make steps towards this goal. 
4.1 Explicit symmetry group of NHDM 

An explicit symmetry of the Higgs lagrangian is a transformation of the doublets <p a that 
leaves the lagrangian invariant. This is to be contrasted with the notion of reparametrization 
invariance of the model mentione in the introduction, which refers to the invariance of physical 
observables but not the parameters of the lagrangian itself. In simple cases such symmetries 
are obvious and can be revealed directly from inspection of the Higgs potential. A more 
intricate situation takes place when symmetries are hidden. They could become obvious after 
an appropriate reparametrization transformation is applied, but such a transformation cannot 
be guessed in any simple way. 

There are several issues concerning possible explicit symmetries of NHDM. In this paper 
we study the following questions: 

• Which explicit symmetry groups are allowed in NHDM? 

• How can one find the explicit symmetry group of a given NHDM potential? 

• How can one construct examples of NHDM realizing a given symmetry? 

Here we make steps towards answering the first two questions, and the answer to the third 
one will become clear on the way. 

We tackle the problem in the adjoint space R , which arguably gives the clearest view of 
possible symmetries of NHDM. Since we want the kinetic terms to be invariant, we consider 
here only (anti)unitary reparametrization transformations, which induce proper rotations of 
r, and keep tq and r 2 unchanged. Therefore, any symmetry of the potential is automatically 
a symmetry of the entire Higgs lagrangian. 

Let us first derive a general expression for the symmetry group of an NHDM potential. 
We rewrite the Higgs potential ffl5l) in the following form 

V = -M r + -A 00 r^ + (-M + L^n + -A^r, , U = A 0i . (19) 

Let us consider this expression in the entire adjoint space and study symmetry properties of 
vectors Mj, Lj, and of the tensor Ay separately. In the discussion below we use short notation 
d = N 2 - 1. 

Let us denote by Gm the subgroup of 0(d) that leaves the vector Mj invariant. If Mj is 
non-zero, this group is given by all rotations and reflections in the (d — l)-dimensional space 
orthogonal to Mf. Gm = 0(d — 1). If Mj = 0, this group coincides with the entire 0(d). 
The same is valid for Gl, the symmetry group of Lj. Obviously, if Mj and Lj are linearly 
independent, the two 0(d — 1) groups are different subgroups of 0(d) and intersect across 
0(d-2). 

Let us denote by G\ the symmetry group of the tensor Ay. Since Ay is real and symmetric, 
it can always be diagonalized by an SO(d) rotation, and the symmetry group Ga is then defined 
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by the spectrum of Ay. If all d eigenvalues are different, then G\ = (Z 2 ) d . If some of the 
eigenvalues coincide, then it is promoted to a continuous group according to the obvious rule: 
each fc-fold degenerate eigenvalue changes (Z 2 ) k —> 0(k). The generic representation is then 

G A = 0(h) x 0(k 2 ) x ■ ■ ■ x 0(k p ) x (Za)"-* 1 --*" . (20) 

Let us now describe the intersection of these three groups, Gm H Gl H Ga, which gives the 
symmetry group of the quadratic form (fT9|) in the entire space R d of vectors r. This symmetry 
group consists of transformations, which leave invariant each of My Li, and Ay. This group 
is non-trivial if and only if there is an eigenvector of Ay orthogonal to both M» and Lj. This 
statement can be also put in algebraic terms. Let us introduce the short notation A n a for a 
vector that can be obtained by acting of the n-th power of tensor Ay on a vector cij. Then, 
let us denote the external product of d vectors dj, fej . . . Zi as 

(o, b, . . . , z) = £iii2---id ' ' ' ) 1 ) 

where eyy...i d is the standard fully-antisymmetric tensor in d dimensions. Finally, let us 
introduce the r -dependent vector fii(r ) = — Mi + Lir . Then, the necessary and sufficient 
condition for existence of an eigenvector of Ay orthogonal to both Mi and Lj is that the 
following polynomial in powers of tq is identically zero: 

X(r ) = 0, Afi, A V, ■ • • , A d ~ V) = • (22) 

If needed, this single criterion can be rewritten in terms of M, and Lj separately. 

If this criterion is satisfied, Gm H Gl fl Ga is a non-trivial group. How large this group is 
depends on the number of eigenvectors of Ay orthogonal to both Mi and and on degeneracy 
of the corresponding eigenvalues. Each such eigenvector contributes a factor of Z 2 to the 
group, and each fc-tuple of such eigenvectors with degenerate eigenvalues gives rise to a factor 
of O(k). The resulting expression has the form of (|20|) . with the total dimension d replaced 
by the dimension of the eigenspace orthogonal to Mj and Lj. 

The last step in identification of the symmetry group of NHDM is to take into account that 
not all orthogonal transformations in 0(d) are realizable as reparametrization transformations 
of the doublets. This group of physically realizable transformations, which we denote by Go, 
coincides with the symmetry group of the orbit space. It consists of the reparametrization 
group SU(N) together with the corresponding anti-unitary transformations in the adjoint 
representation. With this notation, the symmetry group of NHDM can be finally written as 

G = G M n G L n G A n G . (23) 

Thus, the question of identification of all possible explicit symmetries of NHDM is reformulated 
as the search for all possible types of intersections of the individual groups in ( 12"3"|) . We postpone 
the complete analysis for a future study. 

We stress that even if identification of the symmetry group of a given Higgs potential is 
difficult, the inverse task — construction of potentials with desired symmetries — is much 
easier. Using the freedom in choosing M M and A Mi ,, we can align the vectors and the tensor 
with the orbit space in the way we like. The reader will find examples in Section |5j 
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4.2 Which discrete groups are realizable? 

Additional discrete symmetries play an important role in theories beyond the Standard Model. 
However, not all of them are equally interesting in the context of our analysis. Indeed, some 
discrete group will shape the phenomenology of the model in a distinct, group-specific way. 
Other discrete groups, for example, any discrete subgroup that arises as a subgroup of a 
continuous symmetry group, do not tell us about the Higgs potential anything specific to this 
group (see an extensive discussion of this point in [26J). 

The distinction is made clear by the following definition: we call a discrete group H 
realizable (in NHDM) if there exists an NHDM Higgs potential symmetric under H, but not 
symmetric under any continuous group containing H. For example, the NHDM potential with 
cubic symmetry, 

N . N 

V = -m 2 5>l0 a ) + \ ^(0l0 a ) 2 , (24) 

a=l a=l 

contains the symmetry group Sjq of permutations of the doublets, which is not a subgroup 
of any continuous symmetry group of this potential. Thus, is realizable in NHDM. On 
the other hand, the same potential is also invariant under phase shifts of the first doublet by 
multiples of 2ir/q, for any integer q, giving a specific realization of the Z q symmetry group. 
However, this group is not realizable because it arises not on its own, but as a subgroup of 
the full U(l) symmetry group of arbitrary phase shifts of the first doublet. 

A natural question is then: which discrete groups are realizable in NHDM? Note that the 
groups entering ( 123]) contain either continuous groups or factors of Z 2 (which makes Z 2 a 
realizable group). Nevertheless, we have just given an example of a realizable discrete group 
different from ■ Clearly, such symmetry groups must arise from non-trivial intersections 
of continuous groups. 

Let us illustrate this point with a specific example. In 3HDM, the field space is invariant 
under cyclic permutation of the three doublets <pi — > — > 03 — > <f>%. However, this Z 3 group 
is a subgroup of the U(l) group generated by — A2 + A5 — A7. The corresponding finite-angle 
rotations of the doublets, <p a = R a b<Pb, are given by the matrix 

1 / 1 + 2 cos a 1 + 2 cos a" 1 + 2 cos a' \ 
R a = - I 1 + 2 cos a' 1 + 2 cos a 1 + 2 cos a" , 
3 \ 1 + 2 cos a" 1 + 2 cos a' 1 + 2 cos a J 

The three elements of the Z3 subgroup correspond to a = 0, 27r/3 and 47r/3: 

/ 1 \ / 1 \ / 1 \ 

#0= 010 , R.2* =001, i?4.= 100 . (26) 
\ 1 / \ 1 / \ 1 / 

In the adjoint space, the transformation (|25p with a generic a involves all eight coordi- 
nates rj. There are two invariant subspaces of vectors f under the action of this group: 
-E13468 = T3, r4, Tq, Tg) and E257 = (r 2 ,r 5 ,r 7 ). However, at a = 27m/3, the first invariant 
subspace can be split into two smaller invariant subspaces: the root plane E 3 s = (r 3 ,rs) and 
its complement E ue = (ri,r 4 ,r 6 ). In fact, upon the action of so-defined Z 3 group we have 
precisely the cyclic permutations among the "usual", "primed" and "double-primed" sets of 
coordinates discussed in Section El 



2ir Aw 
a' = a + — , a" = a + — . (25) 
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Let us now consider a Higgs potential with a maximal orthogonal symmetry in each of 
these subspaces, E 38 , E U6 and £257, i>e. with symmetry group Gm H Gl H Ga = Ga = 
0(2) x 0(3) x 0(3): 

V = -M r + ^A r 2 + ^A 38 (r 3 2 + r 8 2 ) + ^A 146 (r 2 + r\ + r 2 ) + ^A 257 (r 2 + r 2 + r 7 2 ) , (27) 
with suitable values of M and A's. Here we chose M{ = and 

A-fiv = diag(Ao, A146, A257, A38, Ax46, A257, Aug, A257, A38) . 

Due to A38 7^ A146, this potential is not symmetric under the entire U(l) group (|25l) . but it is 
symmetric under the discrete Z3 subgroup (|26l) . Thus, this subgroup arises from intersection 
of the symmetry group of the orbit space, Go, and the symmetry group of the potential, which 
in this case is given just by Ga. 

With this example, it becomes clear how to identify a realizable discrete group in NHDM 
(other than powers of Z 2 ): we just need to find a discrete group H whose invariant-subspace 
decomposition of the r-space differs from the G-invariant-subspace decomposition for any Lie 
group G C G that contains H. 

Indeed, let us pick up a Lie group G C G , which contains a discrete group H, and let 
the G- invariant-space decomposition of the r-space be @Ei. Clearly, the if -invariant-space 
decomposition, ®Ei )(X , can be only finer, not coarser; that is, it can be obtained form ®Ei 
by either keeping each E{ intact or splitting it further into smaller subspaces Ei — > ® a E^ a . 
Now, recall that we have full freedom when constructing the Higgs potential. In particular, we 
can construct a Higgs potential whose symmetry group is exactly <S>0(ki i(X ), where each factor 
0(ki ya ) is the symmetry group of rotations inside the if-invariant space Ei >a of dimension fcj jQ ,. 
If the H- and G-invariant-space decompositions are not identical, then the potential is not 
symmetric under the full group G. If this situation holds for all possible Lie groups G C Go, 
then the discrete symmetry of the potential under H cannot be extended into any continuous 
one, thus, the discrete group H is realizable. 

This strategy can be applied to NHDM to generate the full list of discrete groups realizable 
at given N. We postpone this analysis for a future work. 

4.3 Spontaneous symmetry breaking 

Even if the Higgs lagrangian has a non-trivial explicit symmetry, this symmetry can be spon- 
taneously broken in the vacuum state. In other words, the symmetry group in the vacuum, 
G v , can be smaller than the explicit symmetry group G of the Higgs potential. 
This possibility brings in several natural questions: 

• How exactly can each of the explicit symmetries be spontaneously broken? 

• How can one recognize at the level of potential that an explicit symmetry is sponta- 
neously broken? 

• How can one construct examples of NHDM realizing a given pattern of spontaneous 
symmetry breaking? 
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Let us first give a general expression for G v . In the derivation of the explicit symmetry group 
G given by Eq. [231 we searched for all possible orthogonal transformations of the orbit space 
that left the potential invariant. Now, deriving G v , we require that only those transformations 
that leave the v.e.v.'s intact are allowed. In other words, we must replace Go by the isotropy 
group of (f), which we denote by G^: 

G v = G M n G L n G A n G (f) . (28) 

Using the isotropy groups for the charge-breaking and neutral vacua found in |29j, one can 
in principle write down the full list of spontaneous symmetry breaking patterns possible in 
NHDM. We postpone this study, too, for a future work. 

If we know that the Higgs potential has a non-trivial explicit symmetry, i.e. (122j) is satisfied, 
and if we have an expression for the position of the global minimum (r M ), then we can check 
whether the vacuum violates this symmetry. Geometrically, the criterion is whether the vector 
(j-j) is orthogonal to the same eigenvectors of as M, and Lj. If this does not hold for at 
least one eigenvector, then the symmetry is spontaneously broken. 

These observations are straightforward, but they do not resolve the problem of detecting 
spontaneous symmetry breaking at the level of potential (that is, they require the knowledge of 
(r M )). In principle, if we know the potential, then all the properties of the vacuum including the 
possible spontaneous symmetry breaking must follow. There must exist a condition involving 
only A M j, and M p for spontaneous breaking of a given symmetry. For example, in 2HDM this 
condition for spontaneous breaking of CP- violation was derived in [151 US] an d was based 
on non-trivial geometric properties of the potential. It remains to be seen whether a similar 
derivation can be conducted in a general NHDM. 

5 Examples of 3HDM with specific symmetry patterns 
5.1 Decoupling of directions 

In this Section we will make use of the general results of the two previous Sections and construct 
several 3HDM potentials, which exhibit various patterns of discrete symmetry breaking. The 
strategy we use is the following. In each case we first consider a lower- dimensional section of the 
orbit space, and then draw a quadric (the critical equipotential surface) with desired symmetry 
that touches the orbit space at isolated points. Then we write down a Higgs potential that 
corresponds to such an equipotential surface; the symmetries of the Higgs potential and of the 
vacuum states will follow. 

Two issues in this procedure require justification. First, let us reiterate the point that 
when constructing the Higgs potential we have the full freedom in choosing its parameters, 
provided the positivity conditions are satisfied. It means that any quadric, no matter where 
it is placed and how it is oriented with respect to the orbit space, is an equipotential surface 
of a suitable Higgs potential. 

The second point is that when constructing specific cases we can safely focus on any affine 
subspace of the adjoint space M N without worrying what happens outside it. In order to 
illustrate this procedure of "decoupling" of unwanted directions, let us choose a diagonal and 
non-singular A^. The linear term in the potential can be compensated by a shift of (ro,rj), 



15 



Figure 4: Illustration of decoupling of unwanted directions, see explanations in the main 
text. The gray circle shows schematically the orbit space, the hyperboles are the critical 
equipotential surfaces. 

and up to the overall constant we have 

V = h (r - c ) 2 + l - Aifa - c,) 2 • (29) 

i 

Here, Co = M Q /X Q , Ci = Mij\. All c's and A's can be chosen independently. Suppose that 
A > 0, and one of A,, for example, Ai, is positive and very large: Ai 3> Ao, |A»|, i ^ 1. 
Then, the equipotential surfaces are very oblate in the direction of n, see Fig. 0J left. An 
equipotential surface M. c with a sufficiently small C approaches the orbit space most closely 
at ri — > c\ and quickly runs away as r± deviates from c\. In the limit Ai — > oo, all the contact 
points between the equipotential surface and the orbit space, which define the minima of the 
potential, are restricted to the subspace r 1 = c\. So, one can focus on this subspace and 
disregard the fact that the orbit space extends beyond it. In this way several directions can 
be "decoupled" from the analysis. 

This procedure can be called the "strong asymptotic decoupling" of unwanted axes. It 
works for any contact — symmetric or asymmetric, from inside or from outside — between an 
equipotential surface and the orbit space. However, in certain cases a relaxed and more useful 
version of decoupling can be used. For example, let us consider a constant-r section of the 
orbit space. The lightcone condition restricts the orbit space to lie in a region on and inside 
the unit ball Yli 71 ! — 1- Let us again pick up axis 1 and consider various n\ = const sections 
of the orbit space. The orbit space has the largest extent along n,, i ^ 1, at ri\ — 0, which is 
very schematically illustrated by Fig. HJ right. Therefore, if c\ = in (129"j) . if it sufficient to 
take any non-negative Ai to safely decouple axis 1 and guarantee that the global minimum is 
located at n\ = 0. This includes the "cylindric case" , Ai = 0, i.e. when the potential does not 
depend on the coordinate n\ at all. 

5.2 Spontaneously broken Z 3 : neutral vacuum 

We now focus on 3HDM and start with the Z% symmetry group generated by cyclic permuta- 
tion of the three doublets, discussed in Section 14.21 A Higgs potential explicitly realizing this 
symmetry is: 

V = -M Q r + ^Aor 2 + iA 38 (r 2 + r 2 ) + ^A ± (|r 12 | 2 + |t- 45 | 2 + |r 67 | 2 ) . (30) 
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In terms of doublets, the potential has the form 



V 



-^[(0 t i0i) + (402) + (^0 3 )] + 



An + A 



■38 



6 



+ (4<t> 2 y + (44> 3 y 



+ 



2A n - A 



^-^[(0t 01 )(0t 02 ) + (0t 02 ) (0 t 03) + (0t 03 ) (0 t 0l)] 



+ Y [l^2| 2 +|0^ 3 | 2 +|030ll 

Note that in principle we could add terms linear in r± + r 4 + r 6 and r 2 + r 5 + r 7 , which are also 
■^-symmetric. However, we prefer to use the geometric argument just discussed to decouple 
all the "transverse" directions and to focus only on the root plane (713, ng). 

To spontaneously break the ^-symmetry, the potential must have the global minimum at 
three points lying on the root plane: 
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The critical equipotential surface is therefore the circumscribed circle of the triangle, see 
Fig. El left. In order for these points to be minima, it is sufficient to take any A^ > 0, together 
with positive A and negative A 38 , provided that A + A 38 > 0. In the entire r M space, the 
critical equipotential surface is a hyperboloid that touches the forward lightcone along the 



circle r\ + r§ 



ri at some tq. 





Figure 5: Spontaneous Z 3 -violation into neutral (left) and charge-breaking (right) minima. 
The critical equipotential surfaces are circles; the points of the global minimum of the potential 
are indicated by dots. 



Expanding the potential at a minimum point up to quadratic terms, one can find the Higgs 
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boson masses: 



m l = yOo + A 38 ) , m 2 H± = — (-A 38 ), rn 2 H0 = — (A± - A 38 ) . (31) 

Here /i is the Higgs boson along the v.e.v., H° are all the remaining neutral Higgs bosons, and 
H ± are all the charged Higgs bosons. 

Note that in addition to the discussed ^-symmetry, this potential has other explicit sym- 
metries, both discrete and continuous. We do not discuss them here because the sole purpose 
of this example was to illustrate the Z 3 breaking in the computationally simplest form. 



5.3 Spontaneously broken Z 3 : charge-breaking vacuum 

Consider now the same Z 3 -symmetric potential ( 130]) . but with A 38 > 0. In this case the 
equipotential surfaces in the r M -space are concentric ellipsoids sitting inside the "hole" of the 
orbit space. The critical equipotential surface is the ellipsoid that touches the orbit space from 
inside. If Aj_ is large enough, then the contact points are again located in the root plane. The 
(n 3 ,n 8 ) section of the critical equipotential surface is represented on this plane by the circle 
of radius 1/2, see Fig. [5] right. 

This potential has the global minima at the three maximally charge-breaking points: 



«" : <^7i(o)' <*>=(!!)• 

An explicit calculation gives the following Higgs masses: 

2 A_i_ 2 Aj_ — A 38 2 A 38 2 4A + A 38 2 . . 

mj = — v , v , — v , v . 32 

2 4 2 12 

The positivity of all the masses squared are guaranteed by Aj_ > A 38 > 0. 



5.4 Other symmetry patterns on the root plane 

Other patterns of symmetry breaking with minima localized on the root plane can be devised 
with Higgs potentials having less symmetry in the root plane than (130]) . This can be achieved, 
for example, by replacing 

^3 8 (r 3 + rl) -> A 3 rJ + A 8 r^ with A 3 ^ A 8 , (33) 

and by introducing M 3 ,M 8 ^ 0. The (n 3 ,n 8 ) sections of the equipotential surfaces can then 
be ellipses, hyperboles or parabolas, depending on the signs of A 3 , A 8 (as usual, we keep A 
and A^ positive). Several examples are shown in Figures [6] and [0 
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Figure 6: Spontaneous ^-violation into neutral (left) and charge-breaking (right) minima. 
Hyperboles are the critical equipotential surfaces, the positions of the minima are indicated 
by dots. 



Figure 7: Coexistence of degenerate minima with different symmetries. Left: coexistence of 
■^-symmetric and Z2- violating minima. Right: coexistence of charge-breaking and neutral 
minima. 

• A 3 < 0, A 8 > 0: spontaneous breaking of the Z 2 symmetry into neutral (Fig. EJ left) or 
charge-breaking (Fig. [61 right) vacuum. 

• A 8 < A 3 < 0: triple degenerate neutral vacuum, Fig. [TJ left, with coexisting points that 
break and conserve the Z2 symmetry. A similar picture for charge-breaking vacua is also 
possible for Ag > A3 > 0. 

• A3 > 0, A§ < 0: double degenerate minimum with coexistence of neutral and charge- 
breaking vacua, Fig. right. In this Figure we deliberately shifted the equipotential 
surface to illustrate that the coexistence and degeneracy of neutral and charge-breaking 
vacua does not require any symmetry of the potential. 

In all cases the values of M3 and Mg are assumed to be appropriately adjusted. 




n 
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5.5 Spontaneous breaking of Z} 4 



More complicated symmetry breaking patterns arise in the transverse space. For example, let 
us construct a Higgs potential with the dihedral symmetry group the symmetry group of 
a square. A particular representation of this symmetry in the space of Higgs doublets is given 
by the group of sign flips of the first or of the second doublets as well as exchange 0i -H- 2 
with 03 always kept fixed. In the adjoint space, these different kinds of Z 2 transformations 
correspond to: 

01 -)■ -0i n X2 -)■ -ni 2 , n 45 -)■ -n 45 , 

02 -> -02 "12 -> -Tli2 , n 67 -»■ -n 67 , 

0i ^ 02 => "3 ->■ -w 3 n 2 -» -ra 2 n 45 <H- ng 7 . (34) 
An example of a ^-symmetric potential is 

V = -Moro + ^Xorl + ^Xxrl + ^X^rl + r^ + ^jirl + rl + r 2 ^ (35) 

= -^[(4<h) + (0 2 2 ) + (0303)] + + (4^2) + (0 3 03)] 2 

+ ^[Re(0l0 2 )] 2 + % |[Re(0 2 3 )] 2 + [Re(0 3 0i)] 1 



2 

-^257 



+ 2 



{[Im(^V 2 )] 2 + [lm(^0 3 )] 2 + [Im(0S0!)] : 



with positive Ao, Ai and A 2 57, and negative A46. Using the decoupling arguments, we can be 
sure that the global minimum is located in the ni = n 2 = = nj = subspace. Then we 
note that the Higgs potential ( 135]) does not depend on 113 and n 8 ; therefore, the equipotential 
surfaces are cylindric along these two directions. Projecting them onto the (n^ne) plane, we 
find concentric circles around the origin because the potential is 0(2)-symmetric in this plane. 

A b A b 





Figure 8: Patterns of symmetry violation on the (714, ne)-plane. Left: spontaneous violation 
of -D4 down to Z 2 ; right: spontaneous violation of Z 2 x Z 2 . 



In the Appendix we describe the geometry of the n\ = n 2 = = n 7 = section of the 
orbit space. When projected onto the (714, n 6 ) plane, it has the square shape shown in Fig. El 



20 



The four vertices correspond to neutral vacua with coordinates 

, Vs Vs 1 

n 4 = ±— , n 6 = 0, n 3 = — — , n 8 = --, 

n 4 = 0, n 6 = ±—, n 3 = +—, n 8 = -- . (36) 

The v.e.v.'s of the lower components of the three doublets at these points are 

<$>,<$>) = ^(0, ±w, v) or -L(±v,0,v). (37) 

The critical equipotential surface of the potential ( 131)]) is the circle that goes through the 
four vertices of the square, Fig. [BJ left. One gets four degenerate global minima, and at each 
minimum the D4 symmetry group is broken to Z2 group of sign flips of the doublet with zero 
v.e.v. Thus, we have an example of spontaneous breaking D4 — > Z2. An explicit expansion of 
the Higgs potential at the minimum gives the following physical Higgs boson masses: 

hf: m? = -^V, -^V, (38) 

,0 2 ^1 2 2A257 — -^46 2 ^46 2 ^Ao + 3A46 2 4A 2 57 — A46 2 

hi- m^—v, v , — — v , v , v . 

All these masses squared are positive once > A46 > — 4Ao/3. 

One can also think of a less symmetric potential in the (714, 7Z6)-plane. For example, by 
replacing 

A 46 (r4 + rl) -> A 4 6+(r 4 + r 6 ) 2 + A 46 _(r 4 - r 6 ) 2 , (39) 

with A 46+ 7^ A 46 _ but both being negative, we get a potential symmetric under the group 
Z2 x Z2 generated by <pi <H- <p2 and by simultaneous sign flips of <fi\ and </> 2 (with kept 
constant). Its equipotential surfaces on the (n 4 , n 6 ) are ellipses elongated along the n 4 + n 6 = 
or n 4 — hq = directions. The critical equipotential surface will again touch the square at 
four vertices as shown in Fig[HJ right. The Z2 x Z2 symmetry is broken completely at each of 
the four degenerate minima. 

5.6 Spontaneous breaking of the octahedral symmetry 

Other symmetry groups can be realized in higher- dimensional subspaces. Let us consider, for 
example, the following potential 

V = -M r + ^\ rl + ^\246(rl + rl + rl) + ^\ 157 (rl + rl + r 2 7 ) (40) 

= - + (4^ 2 ) + (4<Ps)} + y[(40O + (40 2 ) + (4fo)] 2 

+ ^f- {[Im^)] 2 + iMtifo)} 2 + iM4<Pi)f 

+ {[Re(4^)] 2 + [lm(0 2 3 )] 2 + M^i)] J 

with positive A , Ai 57 and negative A 2 46- The mismatch between the real and imaginary 
parts of the doublets' cross products is introduced here on purpose. Due to the decoupling 
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argument, the global minimum is located in the n\ = n 5 = n-j = subspace. Within this 
subspace, the equipotential surfaces are cylinders along 713, n 8 and are concentric spheres in 
the (n 2 , n 4 , n 6 )-subspace. 

In the Appendix we show that the n\ = = nj = section of the orbit space, when 
projected into the (n 2 , n^, n 6 )-space, has the shape of a slightly convex octahedron, which is 
schematically shown in Fig. |9j The vertices and the edges are indeed those of an octahedron, 
while the faces are not planar but slightly convex. 




Figure 9: Spontaneous breaking of the octahedral symmetry in the (n 2 , n^, r^-space. The 
octahedron drawn here depicts the orbit space in this subspace, and it has slightly convex 
faces. The critical equipotential surface not shown here is the large sphere that touches the 
six vertices of the octahedron. 

Within the (n 2 , n 4 , n 6 )-subspace, the octahedral symmetry group Oh arises from permuta- 
tion of the three axes n 2 ,n 4 ,n 6 as well as their sign flips. This group can be extended into 
the entire adjoint space. At the level of Higgs doublets, it is generated by the following trans- 
formations: the CP-transformation (which flips axes n 2 , n^, n-?), the sign flips of individual 



doublets, the 
with matrix 



Z2 transformations such as 2 -H- i03, and the Z3 transformation 



R 



ab 












ab(Pb 



(41) 



which generates the cyclic changes n 2 — > —ne — > — > n 2 and ri\ — > nr — >■ — > n\. Clearly, 
the Higgs potential f l4"0]) is also symmetric under Oh- 

The critical equipotential surface of the potential is the circumscribed sphere of the octa- 
hedron, which touches it only at the six vertices labeled in Fig. [9] as P 2 , P4, P§. These vertices 
correspond to the neutral vacua with the following v.e.v.'s of the lower components of the 
doublets: 



«0?>, <$» 



V2 



(v, ±iv, 0) 



or 



V2 



(v, 0, ±v) 



or 



V2 



(0, v, ±v) 



(42) 



Thus, there are six isolated degenerate minima of the potential ( 1401) . and at each of these 
minima the symmetry group Oh is broken down to P4. An explicit expansion at any of these 
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six minima gives the following physical Higgs boson masses: 



hf: "? = -^fV, ~^fv\ (43) 

,0 2 ^157 2 ^157 — -^246 2 "^246 2 ^Ao + 3A246 2 ~-^246 2 
h { : 7T7,j = — — V , V , —V , V , V . 
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Notice that in ( H2"j) two minima out of six have relative phase between the two non-zero 
vacuum expectation values. However this does not imply that these vacua are CP-violating 
H. Indeed, these minima are physically equivalent to the other four, so that there exists a 
basis transformation from Oh that rotates away the phase without changing altogether the 
parameters of the potential. Thus, despite a relative phase, all six minima are CP-conserving. 
Examples of this situation are also known in 2HDM, see e.g. [28]. Whether 3HDM permits 
coexistence and/or degeneracy of CP-conserving and spontaneously CP-violating minima 
remains to be studied. 

Similarly to the previous case, less symmetric situations can be easily devised by consid- 
ering the ellipsoidal rather than spherical equipotential surfaces. 



6 Discussion and conclusions 

The principal result of this paper is development of a way to efficiently treat the general 
NHDM Higgs potentials. We illustrated this approach with several particular 3HDM exam- 
ples, which exhibited various symmetry breaking patterns. We stress that these examples do 
not exhaust all the possibilities in 3HDM; however, they revealed a number of interesting dif- 
ferences between the scalar sectors of 2HDM and 3HDM. In particular, we made the following 
observations about the minima and symmetries of 3HDM potentials: 

• minima that spontaneously break and conserve a Z2 symmetry can coexist and even be 
degenerate, 

• neutral and charge-breaking minima can coexist (this was also noted in [24J) and be 

degenerate, 

• there can be triple-degenerate charge-breaking vacua and six-fold degenerate neutral 
minima, 

• degenerate minima do not necessarily imply that the potential itself has any symmetry. 

The last statement comes from the simple observation that we can break the ^-symmetry of 
the potential f l3"U|) by drawing an ellipse in the root plane and placing it asymmetrically with 
respect to the triangle, so that it still passes through all three vertices of the triangle. The 
three neutral points will then be the degenerate global minima of the potential, although with 
distinct Higgs mass boson masses. This should be contrasted to 2HDM, where in the tree-level 
approximation degenerate minima can occur only as a result of spontaneous breaking of some 
symmetry, [15J. 

In principle, all the specific 3HDM examples we considered could be easily analyzed with 
straightforward algebra, which we indeed used to calculate the Higgs boson masses. However, 

2 1 am thankful to Pedro Ferreira for this remark. 
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it is through the geometric methods that these examples were found. The same geometric 
arguments also guaranteed that the minima we considered were global. However, the power 
of geometric approach goes beyond such simple examples, as one can easily imagine less 
symmetric types of contact between the critical equipotential surface and the orbit space. 

For example, in the symmetric cases we had several degenerate global minima in the orbit 
space. However this degeneracy can always be lifted by a sufficiently small non-symmetric 
shift of the critical equipotential surface. In this way we can get, for example, a potential 
with the charge-breaking minimum lying slightly higher or lower than the neutral one, or a 
potential with several neutral minima lying at different depths. 

Let us also note that in all the cases we considered the degenerate minima were located 
at different r« but equal ro- This property can be easily broken if a sufficiently skewed (due 
to non-diagonal A^) or displaced (due to non-zero M;) equipotential surfaces are considered. 
However, the experience from 2HDM lends support to the argument that such skewed or 
displaced cases will not generate any new pattern of symmetry breaking in addition to what 
can be generated with more symmetrically placed equipotential surfaces. 

Let us also comment on the number of degenerate minima in 3HDM. Due to the Bezout 
theorem in algebraic geometry, a second-order (the critical equipotential surface) and a third- 
order (the charge-breaking orbit space) surfaces intersect along an algebraic manifold with 
order less or equal to six. If this intersection is meant to be a union of disjoint points, 
one needs to allow two orders per point, giving thus three points. Therefore, the charge- 
breaking global minimum has degeneracy in 3HDM not more than three. Our example with 
spontaneously broken Z% realizes this maximally degenerate symmetry breaking to the charge- 
breaking vacuum. 

The neutral orbit space is represented by an intersection of the same cubic surface and the 
forward lightcone, which is a quadric. Therefore, the neutral orbit itself is a sextic manifold. 
Its intersection with the orbit space has order 12, which can encode up to six separate points 
in the orbit space. Therefore, the neutral order space can have degeneracy up to six. The 
six-fold degenerate neutral vacuum was realized in the example with the octahedral symmetry. 

One can think of several directions for future work. 

Perhaps, the first goal would be to give a complete classification of explicit symmetries of 
the NHDM potential for different N. Especially interesting are the discrete symmetries, whose 
spontaneous breaking does not lead to extra massless scalars. One can study phenomenological 
manifestations of various discrete symmetries and check to what extent these symmetries can 
be incorporated into the Yukawa sector and what symmetries in the fermion masses and mixing 
can result. 

Then, one should also try refining the algebraic description of the Higgs sector in NHDM. 
This includes a further study of properties of the NHDM orbit space, both its local structure 
and topology beyond what was analyzed in [29], a derivation of compact algebraic criteria 
for the positivity constraints on the Higgs potential, a derivation of basis-invariant criteria 
for existence of a hidden symmetry, etc. It would be also very interesting to see whether the 
Minkowski-space formalism of [T5J [TBI E] can be as efficient in NHDM as it was in 2HDM. 

In the case of general 2HDM, some dynamical properties such as the mass matrix of the 
physical Higgs bosons [91 [33] and thermal evolution of the vacuum Higgs configuration [19] 
have also been addressed. It is interesting to see if these calculations can be repeated for the 
general N. 
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Finally, we hope that the general approach of this paper, based on the equipotential sur- 
face technique in the orbit space, can be generalized to other Higgs models, involving singlets, 
triplets etc. 

In summary, continuing our study of the iV-Higgs-doublet model initiated in [29] , we 
developed a method which helped us analyze the properties of NHDM Higgs potentials. We 
worked in the space of Higgs field bilinears, and the key roles in our approach were played 
by certain geometric constructions in this space, such as the equipotential surfaces and the 
orbit space. We also studied explicit symmetries possible in NHDM and their spontaneous 
breaking. We illustrated the general discussion with various particular cases of 3HDM, which 
revealed a wealth of possibilities which were not realizable in 2HDM. The technique suggested 
here can be used to gain more further insight into properties of multi-Higgs-doublet models. 
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A Orbit space in the transverse space: specific 3HDM 
examples 

In Section [5] we considered several examples of symmetries in 3HDM coming from the "trans- 
verse space", that is from the space orthogonal to the root plane. Here we derive the orbit 
space shapes in two specific cases: taking the n\ = n 2 = n 5 = n 7 = section of the orbit 
space and projecting it onto the (n 4 , n 6 ) plane, and taking the ni = n 5 = n 7 = section and 
projecting it into the (n 2 , n^, n 6 ) space. We use these projections (disregarding the values of 
n 3 , n 8 ) because the D4- and 0/,,-symmetric potentials considered in Section |5] do not depend 
on n 3 and n 8 . 

We start with the (i-condition ffTT]) and the lightcone condition 

-2 I ,2 , 1 , ,2 , 1 U |2 , 3(p 2 +^ 2 +K 2 )-l 

n =\n 12 \ +|n 12 | +|n 12 | H <1. (44) 

From the fact the any 2x2 minor of the i^-matrix is non-negative, we derive that 

Sp'p" > \n 12 \ 2 , 3pp' > K' 2 | 2 , 3pp" > \n' 12 \ 2 . (45) 

Let us now restrict ourselves to the subspace n\ = n 2 = = n-j = 0. We are interested in 
the projection of the orbit space onto (714,7^) . The (i-condition then simplifies to 

p'n\ + p"n\ = 3pp'p" . (46) 

Let us first consider the cases when at least one among p, p' , and p" is zero. If p = 0, then 
\ n i2\ 2 = \ n 'i2\ 2 = O5 an d we return to the root plane staying on the upper edge of the triangle. 
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After projecting from the (n 3 , n 4 , n 6 , n 8 ) space to the (n 4 , n 6 ) plane, it corresponds to the point 
at the origin. 

If instead p' = 0, then n'/ 2 = 0, so in the "transverse" space the only non-zero component 
is uq, whose value is 

n\ < 3pp" = 3p(l-p). (47) 

As p changes from zero to one, the value of n\ is limited to 3/4, and this particular point 
corresponds to the neutral vacuum. Thus, on the (n 4 , ne)-plane we get a line segment along 
the n 4 = axis from —a/3/2 to v3/2. A similar segment along the tiq = axis arises in the 
case p" = 0. 

Let us now consider the case when all p, p', and p" are non-zero (i.e. the projection on the 
root plane gives a point strictly inside the triangle). Due to 3p'p" > \ni2\ 2 = 0, we necessarily 
stay on the charge-breaking manifold. Then ( 1461) can be rewritten as 



+ —^ = l. (48) 
3pp 3pp 

Thus, for each point inside the triangle defined by p, p', and p", we get an ellipse in the (n 4 , n^) 
plane with semiaxes 

a\ = 3pp' , a\ = 3pp" . (49) 

Note that a 4 + a\ = 3p(l — p); thus the largest possible ellipse with a given eccentricity is at 
p = 1/2. The semiaxes of this largest ellipse can be parametrized as 

\/~3 \/3 7T 

a 4 = — cos 9 , a-6 = si n > < 9 < — . 

We then get a one-parametric family of the largest ellipses with different values of 9. If the 
angular coordinate is a in each ellipse, then 

n 4 = cos 9 cos a , n 6 = sin 6* sin a . (50) 

Thus, the family of ellipses sweep the square on the (n 4 , 7i6)-plane, bounded by 

i v3 (n . a/3 \/3 

n 4 ± n 6 = — coslfr G 



(51) 



Thus, we recover the square of Fig. |8j the vertices and the edges of the triangle correspond 
to the neutral and charge-breaking points, respectively, while the interior can corresponds to 
the charge-breaking orbit space or, if a point lies on an axis, to the neutral orbit space. 

Consider now another section of the orbit space rii = n 5 = n 7 = 0; we are now interested 
in the projection of the orbit space in the (n 2 , n 4 , n 6 )-space. We selected these axes in such a 
way that Re(ni2n' 12 n / 1 / 2 ) = Re(m2n 4 rze) = 0. Thus, the <i-condition simplifies to 

pn 2 + pn 6 + p n 4 = 3pp p . (52) 

Again, a similar analysis shows that the neutral orbit space corresponds to a "cross" , that is, 
three line segments spanning from —a/3/2 to a/3/2 along each of the three coordinate axes. 
If all p's are non-zero, we can rewrite the (i-condition as 



3p'p" 3pp" 3pp' 
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which, at given p's, parametrizes an ellipsoid in the (n 2 , n 4 , n 6 )-space. Considering all the p's 
inside the triangle gives us a two-parametric family of ellipsoids which sweep a certain region 
(n 2 , n 4 , n 6 )-space. This region has square sections in the planes n 2 = 0, n 4 = 0, or n 6 = 0, and 
it has the symmetry of an octahedron. It is convex but the curvature of the faces is sufficiently 
small, so that this region lies inside a sphere with radius a/3/2, touching it only at the six 
vertices corresponding to the neutral orbit space, see Fig. [9j 
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